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The first passage probability (FPP), of trafficked intracellular particles reaching a displace- 
^ ment L, in a given time t or inverse velocity S — t/L, can be calculated robustly from measured 

particle tracks, and gives a measure of particle movement in which different types of motion, 
e.g. diffusion, ballistic motion, and transient run-rest motion, can readily be distinguished in a 
CO single graph, and compared with mathematical models. The FPP is attractive in that it offers 

a means of reducing the data in the measured tracks, without making assumptions about the 
mechanism of motion: for example, it does not employ smoothing, segementation or arbitrary 
thresholds to discriminate between different types of motion in a particle track. Taking experi- 
mental data from tracked endocytic vesicles, and calculating the FPP, we see how three molec- 
ular treatments affect the trafficking. We show the FPP can quantify complicated movement 
which is neither completely random nor completely deterministic, making it highly applicable 
to trafficked particles in cell biology. 

1 Introduction 

^ The ffist passage probability F{t,L), is the probability density per unit time, that a moving 

particle will require a time t before it exceeds a displacement L, for the first time. The concept 
of FPP has a long history in the study of stochastic processes such as Brownian motion and 
^-H chemical reaction kinetics [5]. However, we report here its utility in analysing the trafficking 

of particles in the cell, whose motion is neither random nor deterministic, but complicated by 
many partially understood effects. 

The FPP can be calculated for both observed particle tracks and theoretical models of 
particle motion, for a range of values of the parameter L, and potentially offers a way of 
^ analysing the particle motion, to characterise relevant velocities, diffusivities, or length and time 

scales. Some advantages of using the FPP are that it is a robust measure of particle motion, 
^ and avoids the use of smoothing and segmentation of particle tracks, or the use of thresholds to 

^ discriminate between the different kinds of motion, such as runs and rests. It therefore avoids 

the artifacts that smoothing and thresholding produce. For instance we may be interested 
in the trafficking of a set of fluorescently-labelled particles which are driven intermittantly by 
molecular motors on microtubules and F-actin, but these may also experience advection in the 
flowing cytoplasm, and Brownian motion. A simple way of analysing such complicated motion 
is to reduce its noise by spatial and temporal filtering, and break it up into different components 
[TJ m 121 [51 H] . These procedures are useful and conceptually straight-forward, but have obvious 
disadvantages in that they all involve "fudge factors" , such as smoothing length scales and time 
scales, or thresholds to distinguish how fast, straight or long-lived a track must be in order 
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to count as a segment of motion we want to measure. These procedures weaken the analysis 
since the measured tracks may be far more comphcated than the ideals we have in mind, and 
it is often unclear how much the final result depends on the arbitrary parameters and implicit 
assumptions. 

The present study is motivated by experimental data from fluorescently-labeled Rab5 GT- 
Pase, which plays a role in the coordination of vesicle trafficking in the endocytic pathway of 
eukaryotic cells 7 . Having previously developed a method of tracking hundreds of vesicles 
moving simultaneously in the cell |10j, we began to explore methods of analysing their motion. 
The tracks are characterised by several qualitative features: vesicles make directed runs at a 
range of velocities; these runs persist for variable times and are interupted by rests of variable 
times; after a rest, the vesicles sometimes carry on with motion in the same direction, but often 
make abrupt changes of direction, or seem to reverse their motion along a linear track; and runs 
are sometimes straight or gently curved, but often follow a more erratic path. 

In this article, we demonstrate the calculation of the FPP from experimental data taken 
from control cells as well as cells treated with nozodazolc, which disrupts the microtubules, thus 
compromising the motility of Rab5 GTPase. We compare the results with simple theoretical 
models of intracellular trafficking, and a simple segmentation analysis of the tracks. In a separate 
article, now in preparation, we will present new results on the intracellular trafficking of vesicles, 
using the first passage probability distribution, as well as other analyses. 

2 Experimental section 

2.1 Cell culture and live cell imaging 

HeLaM cells, in DMEM containing 10% FBS, were grown in 35mm glass bottomed dishes 
(MatTek Corporation, Ashland, MA, USA) and imaged at 37°C on an Olympus 1X81 micro- 
scope using a lOOx objective lens, 1.35 NA and an additional 1.5x lens magnifier, fitted with 
a Optoscan high speed dynamic bandpass control monochromator (1800 g mm~^ Holographic 
grating) (Cairn Research, Faversham, UK), and a Photometries Cascade 512 back-illuminated 
camera (Photometries, Tuscon, AZ, USA). Cells were imaged for 500 frames using continuous 
illumination at 10 frames/s. Control experiments showed that imaging under these conditions 
for as long as 50 minutes did not affect the probability that the cells would successfully divide, 
or undergo apoptosis, in the following 24 hours (data not shown). 

Cells were treated with nocodazole, to depolymerise microtubules, as follows. A dish of 
HeLaM cells were incubated in DMEM containing l-/iM nocodazole, at 4°C for 5 minutes, 
followed by 2 hours at 37° C. The cells were then imaged as above. 

2.2 Particle tracking 

Our recently developed tracking method, the polynomial-fit, Gaussian- weight algorithm (PFGW) 
[10', allows us to accurately track the extrema of intensity corresponding to individual vesicles, 
without errors due to the presence of their neighbours in the image. Thresholds of acceptable 
eccentricity, radius, skewness and particle lifetime were employed to reject intensity extrema 
that did not correspond to single vesicles. Manual examination of the tracked movies showed 
that all particles were tracked except for a small minority of faint fast-moving particles. Fig.|2ja) 
shows example tracks for the control cells. 

Static errors may be estimated by calculating the mean square displacement (MSD), as a 
function of time scale t, for the measured tracks [TOl [11]. Fig. [l] shows the MSD for the control 
cells. At the shortest time scale of 0.1 s, there is no plateau that would be characteristic of static 
error. Therefore, the average static error in measuringdisplacements is significantly less than 
•^MSD(0.1 s) = 0.08 fim. It is also apparent from Fig.lllthat the MSD has an exponent greater 
than 1 at i < 1 s and less than 1 at f > 1 s. Thus the vesicle motion is superdiffusive at short 
time scales — i.e. dominated by active motion, and subdiffusive at long time scales — presumably 
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Figure 1: Mean square displacement as a function of time, MSD(t), for all vesicle tracks in the 
control cells. The static error is negligible, and the average motion is superdiffusive at t < 1 s, and 
subdiffusive at t > 1 s. 



because the vesicles tend to stop moving or lose their initial direction after times of the order 
of ~ 1 s. 

2.3 Segmentation analysis 

Apart from the first passage probability calculation, we analysed the tracks by a method of 
segmenting each particle track into discrete runs and rests. This heuristic method is similar 
to methods used to characterise intracellular motion by several previous authors. We briefly 
detail our particular segmentation analysis here. The analysis is based on the observation that 
the vesicles tend to make directed runs along straight or gently curved trajectories, which are 
punctuated by stationary rests. The length, duration and average speed of each run or rest is 
measured by: 

1. A smoothed contour is created to correspond to each particle track. The first step was to 
replace the coordinates of each point in the particle track with the mean of all coordinates 
in the track which were within a threshold distance of Lpi^ = 4 to that point. Thus 
the track is smoothed spatially rather than temporally: this spatial averaging has the 
advantage that short runs of directed motion within a track are more preserved than if we 
took a running average in time of points along the track. The second step was to produce 
a coarse version of the track by successively taking sets of points in the smoothed track 
separated by the distance Lpi^, starting from the initial point, and replacing each set with 
its mean position. Finally, the initial and final displacements in the coarsened track were 
extrapolated so that they extend further than the initial and final points of the original 
track. Examples of track contours can be seen in Fig. |2ja): they consist of straight-line 
sections of length ~ Lpi^. 

2. The position x{T) of a particle along its contour, was calculated as the projection of the 
smoothed track along the contour, taking the projection of its initial coordinates as x = 0. 

3. The parametrised distance x{T) is segmented into discrete runs and rests. Wherever the 
particle moves less than 1 pixel in more than 5 frames, those positions are marked as a 
rest. The remaining sections are counted as runs if the total displacement of each is more 
than 1 pixel. Fig. ^h) shows examples of x{T) for the control data, segmented into solid 
(blue) sections for runs and dotted (red) for rests. 
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Figure 2: (a) Example tracks (white) and smoothed contours (coloured online) of all tracks of total 
contour length greater than 3 fj,m, in one video of the control cells, (b) Examples of parametrised 
distance along the smoothed contours x(T) for the control cells, segmented into runs (solid, blue) 
and rests (dotted, red). 



This procedure ignores other effects observed in the motion of the particles, such as reversals of 
direction and slower erratic motion, and has advantages and disadvantages as mentioned below 
and in our forthcoming article [^|. 



3 Results 

Rab5-labelled vesicles were tracked in control cells, as well as cells subjected to nocodazole 
treatment, as described above. For each case of control or treated cells, at least f videos, of 500 
frames each, were captured of separate cells. The average numbers of particles simultaneously 
tracked at any instant in the control and nocodazole treatment, were 128 and 54 respectively. 

From a set of observed tracks {R„(T)}, the first passage probability distribution is calculated 
by finding the smallest non- negative t that satisfies |R„(T + t) — R„(T)| = L, at each starting 
time point T, at each track in the set, indexed by n. To calculate F{t, L), these values of t are 
counted in bins; then the count is normalised and plotted as a histogram. 

The first passage probability can also be calculated in terms of an inverse velocity: F{S, L), 
where S — t/L is the inverse velocity or "slowness". Since F{S,L)ds — F{t,L)dt, the distribu- 
tions can be converted by scaling t with the length scale: 

FiS,L) = LF{t,L)\t=sL. (1) 

Thus to calculate F{S,L), the set of first passage times t is replaced with S = t/L, be- 
fore binning and counting. Fig. [s] shows F{S,L) for each cell treatment. The distribution is 
normalised by the total number of recorded time points: thus the area under each curve is 
the total probability that a vesicle will make a passage of length L from any starting point. 
It is useful to renormalise this distribution to form a new distribution, Fr{S,L), in which the 
area under each curve is unity. The difference between F and Fr reflects the fact that the 
tracks are finite and often much shorter than the length scales in which we are interested: the 
renormalised distributions only count starting points that result in a passage. The graphs are 
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Figure 3: The first passage probability distribution for all data sets, at a range of L. The distribu- 
tions are normalised (a~b) by the number of data points in all tracks, and (c-d) renormalised by 
the number of passages at each L. The distributions for all treatments show a single peak which 
becomes narrower as L increases. 

qualitatively similar in that they show singly peaked distributions, with maxima in the range 
5=1-10 s /im^^. As L increases, the distributions become narrower, especially in the case of 
the nocodazole treatment. 

We can also plot the first passage probabilities to compare the different treatments. Fig. |4] 
shows F{S^ L) at L =0.5, 1 and 3 /im respectively. As the length scale increases, the difference 
between the treatments becomes more apparent. For example, if we consider an inverse velocity 
of 5 = Is /im~^, corresponding to the peak FPP of the control cells, the probability of 
vesicles traversing various distances, at the corresponding speed of 1/5 = 1 /im s^^, is markedly 
reduced under the nocodazole treatment compared to the control cells (TablejlJ. The differences 
between values of i^(l s /im~^, L) are greater than between Fr{l s /im^^, L), since the proportion 
of vesicles making passages of 1-3 /im is significantly reduced by the nocodazole treatment. 
F{1 s /im~^, 1 /im) is a factor of 10-50 lower for the treated cells than the control cells. In all 
graphs of Fig.[4j we can see that the nocodazole treatment impairs motion at lower S — i.e. the 
probability of particles moving at larger velocities is markedly reduced. The difference increases 
with increasing L: i.e. motion is increasingly impaired at larger length scales. 

In order to test the reproducibility of behaviour between cells, we plot F{S, 1 /im) for the 
tracks from each video from the control data, as well as the combined distribution in Fig. [5j 
The differences in FPP between the control cells are due to the variability in behaviour of cells 
as well as the statistical error of sampling the vesicle tracks. However the differences are much 
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Figure 4: The first passage probability distribution for all data sets, at L = 0.5, 1, 3 fim. Both 
normalisations are shown: F(S,L) and Fr{S,L) (a-c and d-f). Probabilities are markedly reduced 
under the nocodazole treatment compared to the control cells. 
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Table 1: Reference values of first passage probability for all cell treatments. 
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Figure 5: F{S, 1 ^m) for the control cells: from each video (thin lines, coloured online) and from 
all videos combined (solid line H-). 



smaller than those caused by the nocodazole treatment as shown above. 



4 Theoretical models 

The observed first passage probability distributions can be compared with theoretical models of 
various types of particle motion. We consider three simple models below, from which analytical 
expressions for the FPP can be calculated, then expand the discussion to more complicated 
types of motion, in comparison with the experimental data. 



4.1 Steady directional motion 

We consider first the simplest possible kind of particle motion. If all particles were to move at a 
constant speed v, with each moving in any constant direction, then the first passage probability 
would be a delta function, since the time for a particle to move L would always be L/v: 

F{S,L) = S{S-l/v). (2) 

Similarly, if each particle moved at a constant speed in any direction, where its speed was chosen 
from a probability distribution P{v), then 

= ^P{l/S). (3) 

In either case, the first passage probability, as a function of S, is independent of L. This does 
not match F(S, L) from the observed data which becomes narrower with increasing L. Indeed 
steady directional motion obviously does not match the captured videos, in which the vesicles 
show saltatory and erratic motion. 



F{S,L)^P{v) 



dv 
dS 



4.2 Random walks 

For a random walk characterised by a diffusivity D, the probability C(t, r) of finding the particle 
at r at time t is governed by the diffusion equation: 

^^DV^C. (4) 



7 



If the particle is at the origin at t = 0, i.e. C(0, r) = 6{r), then Eq. g can be solved to yield the 
occupational probability C(o)(i,r) at subsequent times [9J: 



(47rD<)»/2 



(5) 



where n is the dimensionality of the random walk. However, we must note that the first passage 
probability is different from the occupational probability: the latter includes the possibility 
that the particle, after it has reached a displacement L, could return to that L. The first 
passage probability can therefore be calculated from Eq. |4] using the same initial condition, but 
imposing a boundary condition C(|R| = L) = 0. This represents a sink of probability wherever 
the particle reaches a displacement of L; C now includes only the trajectories that have not yet 
passed L. In one dimension, these considerations lead to the series solution for the probability 
of finding the particle before the first passage C(bfp) 



(BFP 



COS 



fc=0 



(2A: + l)7rx 
2L 



exp 



(2/fc + l)7r 
2L 



Dt 



The first passage probability is then the probability flux across the boundary: 

F{t,L)^~D(f (VC)-dA, 

J\r\=L 



or in one dimension: 



F{t,L) = 2D 



dC 



dx 



For the one dimensional random walk, the FPP is therefore: 



(6) 



(7) 



(8) 



J^dDRW) {t, L) = jY E(-l)'(2fc + 1) exp 

fc=0 



or in terms of S: 



i^(iDRW)(5, i) = ^ 5](-l)'^(2fc + 1) exp 



fc=0 



(2fc+ l)7r 
2L 



/ (2fc + l)7r 



Dt 



DS/L 



(9) 



(10) 



In Fig. [6] we plot -F'(idrw)('5'i L) for a range of L, taking D — 1. As shown, -F'(idrw)('S'7 L) has a 
single peak which scales as S ^ L. The scaling is evident from Eq. 10 in which i^(iDRW)('S'j L) is 
a function of DS/L. For random walks in two or more dimensions, the first passage probability 
is qualitatively similar, and also scales as 5' ~ L. 

By glancing at Figs. |6]and|3] we can see that the FPP from the measured data bears little 
resemblence to the random walk model, because the peak position in the measured data stays 
approximately constant although L is increased by a factor of 50. 



4.3 Unidirectional runs and rests 

Suppose a particle alternately runs and rests. All runs share a constant direction and velocity 
V, and while running the particle has constant probability a, per unit time, of switching to a 
rest. During a rest, the particle does not move, but has constant probability 6, per unit time 
of switching to another run. These single-particle switching rates a and h lead to a dynamic 
equilibrium where, over a long time, the particle will spend a proportion / = (a/6+ of the 
time running, and 1 — / = {h/a + resting. 

The first passage probability distribution of this "run-rest" model, F(rr)(S', L), should have 
two limits which we can predict before we analyse the model in detail. Over a long track length 
L, the particle spends a total time L/v running, thus the mean number of rests in this period 
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Figure 6: The first passage probability for the one-dimensional random walk, as a function of 
inverse speed, at a range of L, taking D = 1. The peak position scales linearly with L. 

is X = aL /v. If A ^ 1, then, by definition, the first passage probability is dominated by the 
runs: i.e. F(ri^)(S', L) « S{S — v~^) for L ^ v/a. However, if A » 1, the particle experiences 
many rests and runs; as A — *■ cxi, the proportion of time it spends running tends to /. Therefore 
P{i{R){S, L) « 8{S — {vf)~^) for L ^ v/a. Thus we expect to see the first passage probability 
switching from a sharp peak at small-L to another sharp peak at large-L. 

To calculate -F'(rr) in detail, consider a passage of length L that contains n rests. Its total 
time will be 

L " 

i=-+Vrfc, (11) 
V ^ — ' 

k 

if each rest persists for a time . Since there is a constant probability per unit time of the rest 
ending, the probability of that rest surviving for Tk is distributed exponentially: 

Prost(rfe) = be-'^^ , (12) 



and therefore, the probability of a particular combination of rest times {rfc} that satisfy Eq. 11 
is the compound probability 



P({rfc}) = Yl be-^^" = b^'e-''^'-^/^^ , (13) 



k=l 



making use of Eq. [TT] 

The probability of having n rests in the length L follows the Poisson distribution, since they 
occur independently of each other: 



P{n) = , (14) 



where A is given above. 

The probability that the passage requires a time t, given that it contains n rests, is the 



integral of P{{Tk}) over all the possible combinations of {t^} that satisfy Eq. 11 



This integral becomes an integral over a constant once we have substituted Eq. |13[ 

[n - 1)1 
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The probability that the passage requires a time t is then obtained by a sum of P(t\n) over all 
the possible values of n: 

Pit) = E mn)Pin) = e-^Sit - L/v) + ^ ^ ^ ^ lu -lV ' (''^ 

This series solution can be rearranged using the form of the first modified Bessel function of 
the first kind, Ii. 



P{t) = e-^^^^6{t - L/v) + ^^^^e-'*-^^-'^^/^h y2^<^b^ " ^) j " ^^^^ 

Finally, to obtain the first passage probability, we must consider the starting point of the 
particle, since it may be initially running or resting. In the latter case we must include the time 
of the initial rest: 



l-t-L/v 

{t,L) = fP{t) + {l- f) P,,,t{r)P{t-T)dT, (19) 

Jo 

(20) 



^(RR) 



e--v-^(t -L/.) + --" :;;;r^" 

X ( all 



o{2^ab{t-L/v)L/v) + ,^^I,{2^ab{t - L/v)L/v)) , 



or in terms of S: 



(^aIo{2^abL^v-^S - v'^)) + ^ f^h{2,/abL^v-^S - «-i))) 



(21) 



We compare the run-rest model with the observed data by taking the relevant parameters 
of Eq. [2T]from a segmentation analysis of the observed data, as described above. Fig. [7] shows 
statistics of the lengths, times and speeds of runs and rests in the control cells, from which we 
measure v = 1.0 /im s^^, the mean speed of the segmented runs, = 0.76 s, the mean lifetime 
of the segmented runs, and b~^ = 3.8 s, the mean lifetime of the segmented rests. Taking these 
values of v, a and b, i^(RR)(S', L) is plotted in Fig.jS^a). The predicted FPP does indeed switch 
from a delta peak at 5 = ii"^ to a delta peak at 5* = {vf)~^ = 6.0 s /im~^. Having seen 
that Fig. iT^a) shows a distribution of run speeds P{v), we then convolve -^'(rr,) with P{v) to 
compare directly our run-rest model with Fr{S,L) from the measured control data (Fig.jsjb)). 
The run-rest model agrees qualitatively with the measured data, showing a similar narrowing 
with increasing L. However, there is a quantitative discrepancy: the predicted inverse velocities 
are about a factor of 5 lower than the measured values. 

Since calculation of the first passage probability from measured particle tracks is a straight- 
forward procedure, which employs no assumptions about the particle motion, nor arbitrary 
parameters, this comparison between the first passage probability analysis and the run-rest 
model shows that the latter does not accurately refiect the motion of the endocytic vesicles. 
Nevertheless, we do not wish to devalue the segmentation of tracks into runs and rests, because, 
as others have found, it still a useful analysis of particle motion, as it is conceptually simple 
and motivated by direct observation. However, we discuss below some possible reasons why the 
run-rest model is inadequate for predicting the first passage probability distribution of endocytic 
vesicle motion. 



5 Discussion 

We believe that the concept of first passage probability, of particles to be transported a certain 
distance, is a natural way of analysing the trajectories of trafficked intracellular particles, es- 
pecially when our theoretical knowledge of their motion is incomplete. The calculation of FPP 
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Figure 7: Statistics of segmented runs and rests from the control data: (a) run speed, (b) run 
length, (c) run time, (d) rest time. Distributions from of individual videos are shown (thin lines, 
coloured online) as well as the distributions of all videos combined (solid line H- ) . 




model, taking parameters v, a and b 
is evaluated with v = 1.0 /xm s~^. 
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Figure 8: The first passage probability for the "run-rest" 
from the segmentation analysis of the control data, (a) Eq. 

= 0.76 s, and = 3.8 s: -F{rr) switches from a delta function ( — ) at 1/f = 1 s ^m~^ when 
L is small, towards a delta function at l/{vf) = 6.0 s /um~^ when L is large. Colours correspond 
to the legend of Fig. [sl (b) Eq. 21 is convolved with the distribution of run speeds in Fig. 
above, again taking a = 0.76 s, and b^^ = 3.8 s. The predicted F{S,L agrees qualitatively with 
the measurement, showing a narrowing of the distribution with increasing L. 
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reduces the observed tracks without employing assumptions about the nature of each particle's 
motion, nor arbitrary parameters and thresholds. It also leads to a useful graphical interpre- 
tation of measured data: the distribution that results are seen to be qualitatively different for 
conceptually different types of motion, such as random walks, directed motion, and the run-rest 
model. In comparing measured data, the FPP has allowed us to visualise differences in motion, 
in cells treated with nocodazole compared to control cells (Fig. |4] Table [TJ . By comparing the 
FPP, we have seen that each cell treatment has markedly reduced the probability of vesicles 
being trafficked at higher speeds and larger distances. The FPP has also allowed us to visualise 
the reproducibility of motion between different cells (Fig.|5|. 

The FPP can be calculated as a function of different variables. From the point of view 
of the theory, the most natural is time, because for each position R„(T) along a track n, we 
measure the time that the particle takes to make its first passage of length L from that position. 
However, from the point of view of measuring trafficking, in which particle trajectories bear some 
resemblance to ballistic, unidirectional motion, it is natural to plot the FPP as a function of 
inverse velocity S — t/L, or velocity L/t, since either allows us to instantly recognise steady 
unidirectional motion. We choose F{S^ L) for two reasons: inverse velocity is closer to the 
calculation because it only involves scaling the times by L; and a function of inverse velocity 
shows faster moving particles more prominently than a function of velocity, because of the 
implicit transformation between the two (Eq. [ij. 

Three very simple models of particle motion have been compared with the measured FPP. 
The two extremes of steady directional motion and diffusion led to predictions that were qual- 
itatively different from the measured FPP, so a quantitative comparison was unnecessary. The 
run-rest model, of constant run velocity v, and constant "on" and "off" rates, b and a, pre- 
dicted a FPP that was qualitatively similar to the measured data, when we selected a and h 
from a segmentation analysis of the measured data, and allowed a convolution of v over the 
distribution of run speeds from the segmentation analysis. However, the predicted values of S 
were approximately a factor of 5 from the measured distribution. There are, of course, many 
sources of discrepancy between the run-rest model and the experiment; some were mentioned 
above: abrupt changes of direction, reversals of motion along a linear track, curvature of the 
path and slower erratic motion. It is outside the scope of this study to model these effects prop- 
erly: obviously, to model the observed vesicle motion properly would require a much deeper 
understanding of the governing mechanisms. 

However, we suspect that there may be simple interpretations of the coarse discrepancy 
that we have found. For instance, we believe that is wrong to assume there are only two rates 
a and b that govern switching from run to rest and vice versa. For example, a run could be 
interrupted by various different physical processes such as detachment of the motor, obstruction 
of the motor or load, reaching the end of the microtubule, binding of the motor to a regulatory 
protein, etc. We see a few fast-moving tracks that are interrupted by frequent short rests, as 
well as particles which spend most of their time sessile, but make occasional jumps. The fact 
that there must be more than one type of rest is reflected in the non-exponential rest time 
distribution of Fig. [7]^d). 

We must also be aware that it is difficult to assess the reliability of the segmentation analysis, 
especially where the particle tracks do not follow straight lines, and are complicated by several 
types of motion as well as measurement error. For example, our segmentation analysis ignores 
runs that are short in time or length, or slow in speed: these runs would be counted instead as 
part of a rest. Rests that are short in time are also ignored because of the thresholds. Likewise, 
smoothing of particle tracks "irons out" short runs and short rests, making runs and rests appear 
longer in time, and runs appear slower in speed. The process of thresholding and smoothing 
can therefore artificially change our predictions of a, b and P{v) to a significant degree that 
could easily be missed. On the other hand, the FPP analysis is transparent: no smoothing or 
thresholding are involved, and the measured tracks yield a first passage probability distribution 
in a single mathematical step. 
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6 Conclusion 



The first passage probability distribution is a robust measure of intracellular particle tracks 
that allows us to reduce the measured data into a form in which qualitative and quantitative 
differences in motion can be easily visualised and compared. The FPP is transparent in that 
it does not rely on arbitrary parameters or thresholds, nor any implied assumptions about the 
mechanism of a particle's motion, unlike the alternative segmentation analyses. Therefore it 
is an attractive way of analysing intracellular particle trafficking which is neither completely 
random nor completely deterministic. The FPP gives us an illuminating way of comparing 
theories of trafficking with experiment. 
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